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Abstract. We examine the effects of ionization, radiation pressure and main sequence 
winds from massive stars on self-gravitating, clumpy molecular clouds, thereby mod- 
eling the formation and pre-supernova feedback of massive star clusters. We find the 
process of "turbulent mass loading" is effective in confining HII regions. Extrapolating 
typical Galactic high-mass star forming regions to greater initial gas cloud masses and 
assuming steady star formation rates, we determine the timescales for cloud disruption. 
We find that a dense (n c ~ 2 x 10 5 cm -3 ) cloud with initial mass M c ~ 4 x 10 s M 
is able to form ~ 2 x 10 5 Mq of stars (50% efficiency) before feedback disperses the 
gas after ~ 3 Myr. This mass and age are typical of young, optically visible super star 
clusters (SSCs). The high efficiency permits the creation of a bound stellar system. 



1 Introduction 

Most Galactic stars are born in highly clustered regions Jl^|| , where the disrup- 
tive effects of massive stars are paramount. The multitude of dusty high redshift 
sources and the intensity of the far infrared background they produce also imply 
that a major fraction (~ 1/2) of total cosmic star formation has occurred in star- 
bursts, replete with massive stars p0| |. A significant fraction of star formation in 
local starbursts occurs via the creation of super star clusters (SSCs) [ p2p4| , each 
with hundreds to thousands of OB stars crammed into a few parsecs. At least 
some SSCs are gravitationally bound |8|,|| and their masses (~ 10 5 — 10 6 M© 
[ ^9^15( 1) and sizes suggest we may have found globular clusters in their infancy. 

Massive stars violently disrupt their surroundings with ionizing and non- 
ionizing photons, protostellar, main sequence and post-main sequence winds, 
and supernovae. Adding to this complexity is the extremely dense, clumpy and 
turbulent nature of the gas in which high-mass stars are born. We present a 
simplified model of pre-supernova feedback, to examine how the efficient star 
formation required to produce bound clusters may occur in the presence of vig- 
orous energy injection from many massive stars. 



2 Initial Conditions — a Clumpy Molecular Cloud 

We consider a spherical cloud of radius, R c , mass M c , and mean H density n c . 
The cloud consists of dense clumps embedded in a uniform inter-clump medium. 
The clumps are distributed uniformly within a central core of radius 0.2 i? c , 
and with a r _1 distribution outside, which mimics observed molecular cloud 
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profiles. We choose mean clump mass and volume fractions of f m = 0.8 and 
fv = 0.03 respectively. Our clumps, of uniform density n c \, have a mass spectrum 
d A/"d/d In m c i oc ni~^' 6 between upper, to c i, u , and lower, m c y, limits. Clump 
masses span a range of 10 3 , with m c i iU = 0.025A/ C , so no one clump dominates 
the cloud. The clump velocities are set so the cloud is in virial equilibrium. 

Plume et al. [l8j| determined sizes and masses for 25 regions of Galactic high- 
mass star formation. The mean properties of this sample were R c ~ 0.5 pc, 
M c ~ 3800 M Q (virial mass) and thus n c ~2x 10 5 cur 3 (~ 3.2 x 10 4 M Q pc" 3 ). 
Note that these regions are only a small fraction of the host Giant Molecular 
Cloud. Our adopted clump mass and volume fractions imply clump densities 
of n c \ = 6.5 x 10 6 cm -3 and an inter-clump density of 5.0 x 10 4 cm~ 3 . These 
properties form the basis of our model A. Determining the gas properties of the 
precursors to SSCs is more difficult as current observations only probe scales 
down to ~ 30 — 100 pc |?],^3| , while the star clusters have a typical radius ~ 4 pc 
[ p2| . For simplicity we consider models with 10 (B) and 100 (C) times the mass 
of A, but with the same mean density. With the same values of /y and f m 
the clump and interclump densities are also the same. For each cloud we shall 
consider a fixed star formation rate, 050, which converts 50% of the initial cloud 
into stars over 3 Myr - i.e., the approximate time before the first supernova is 
expected to occur. R The properties of our initial gas clouds are listed in table 111 



Table 1. Model Parameters 



Model M C /(M ) fl c /(pc) am/fkms- 1 ) <p 50 /(M Q yr" 1 ) 



A 


4 x 10 3 


0.5 


2.8 


5.8 x 10 -4 


B 


4 x 10 4 


1.1 


5.9 


6.6 x 10~ 3 


C 


4 x 10 5 


2.3 


12.9 


6.7 x 10" 2 



3 Feedback Processes 

To model the feedback processes, we need the total number of hydrogen ionizing 
(A < 912 A) photons emitted per second (S), the bolometric luminosity (Xboi) 
and the mass flux (M w ) and velocity (v w ) of the stellar wind from the forming 
star cluster. We derive these quantities from the STARBURST99 model [|13| 
at solar metallicity. As we are interested only in the first few Myr of stellar 
evolution, we approximate the values of these quantities to be constant for a 
given stellar mass. Motivated by observations of the stellar population of the 
R136 cluster in 30 Doradus 0], we consider a stellar initial mass function (IMF) 

1 This rate excludes the formation of an initial star cluster of ~ 250 Mq necessary to 
produce an ionizing luminosity of 10 49 ionizing photons per second (see below). 
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represented by a broken power law, dA/"*/dlnm* oc (m^/m^o) - ", with a — 1.35 
(Salpeter) for m* > rra*,o = 3 M and a — 0.8 for m* < m*^- We choose 
lower and upper limits of stellar mass of 0.1 M Q and 100 M©, respectively. This 
increases the mass per unit of feedback (assuming negligible contribution from 
stars below m*^) by a factor of 1.35 compared to the standard STARBURST99 
IMF (a = 1.35'between 1 M and 100 M Q ). We find 

S49 = 3.94 x 10 3 M*, 6 , (1) 

ibol,5 = 1-10 x 10 4 Af*, 6 , (2) 

M w _ 6 = 3.09 x 10 3 M*, 6 , (3) 

M w ,-6V 2 ooo = 3.98 x 10 3 M^ 6 , (4) 

where 5 49 = 5/10 49 photons s" 1 , L bo i, 5 = L bo i/1O 5 L , M w _ 6 = M W /10- 6 M Q yr" 1 , 
^,2000 = ^to/2000 km s" 1 and M* )6 = M»/10 6 M Q . The number of m* > 8 M Q 
stars (i.e. number of core collapse supernovae) is A/sn = 1-4 x 10 4 M*fi. 



3.1 Ionization 

Assuming spherical symmetry with the massive stars forming at the center of the 
cloud, we calculate the ionizing flux received at a given distance R from the star 
cluster by accounting for attenuation by clump shadowing, H recombinations 
and dust absorption, so that 

37; = -SA c iAf c i - 47ri? 2 a 2 n i 7i - a d n,iS, (5) 
oil 

where A c \ is the mean cross sectional area of clumps at R, Af c \ is the space number 
density of clumps, a 2 is the recombination coefficient, is the hydrogen nuclei 
number density of the ionized gas and ad = 0.5 x 10 _21 cm 2 0] is the dust absorp- 
tion cross-section per H nucleus. We assume He is single ionized and nuc = O.lnj. 
An HII region rapidly forms with typical size i?st — 3.05 x l0~ 2 Slg 3 n^^ 3 pc, 
for dust free, uniform density gas. The short sound crossing time justifies our 
assumption of uniform rii. Thermal balance maintains a constant ionized gas 
temperature, T,; ~ 10 4 K. Over-pressurized compared to the surrounding neu- 
tral medium, the HII region tends to expand.^] However, when neutral clumps 
become exposed to ionizing photons, they implode and inject mass into the HII 
region [||. A compressed neutral globule remains, which continues to evaporate 
more slowly ||]. We employ models for magnetically supported clumps so that 
the pressure of the ionized gas streaming from the surface of an imploded clump 
at distance R is 

1-8-TxlO^ 1 (^) 4/7 - cl 4/2lKc ^ 3 

= 9.9 x lO^ 21 (^) 4/7 ra ^ 1 Kc m - 3 , (6) 

2 For the propagation of ionization fronts into the clump and interclump material, we 
assume the neutral gas has a temperature of 80 K. 
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where m c \ is the clump mass in solar mass units and n c \ is the initial clump 
density. Clumps may only photoevaporate if this pressure is greater than the 
wind ram pressure (below). Evaporating clumps are gradually ejected from the 
HII region by the rocket effect. However, clumps orbiting in the cloud potential 
continue to enter the HII region and inject fresh gas. The HII region is limited 
by a recombination front as the ionized gas flows out at approximately its sound 
speed. We have examined the details of this mechanism of "turbulent mass load- 
ing" and its implications for Galactic ultracompact HII regions elsewhere f2lf| . 



3.2 Winds 

Stellar winds carve out a hot, low density cavity at the center of the HII region. 
The wind's thermal pressure is suppressed by cooling caused by mass injected by 
photoevaporating clumps (l7). Instead, we calculate the edge of the wind cavity 
by balancing the thermal pressure of the HII region with the wind ram pressure, 

= 7 . 6 x i 6 M -g20oo K cnr3 „ 3 x 1Q 9^ K (7) 

In addition, this pressure contributes to the clump dynamics and may quench 
photoevaporation, particularly at small distances from very massive clusters. We 
treat photoevaporation mass injection into the wind via a supersonic mass-loaded 
wind model, e.g. pj. In cases of extreme mass loading the ionizing photons may 
be trapped in the wind. 



3.3 Radiation Pressure 



Our fiducial clumps absorb most of the momentum flux of the star cluster's 
radiation. The pressure is 

=8 .7x 10 9 ^Kcm- 3 . (8) 
k R pc 



4 Numerical Method for Modeling SSC Formation 

With an N-body code, we follow the dynamics, masses and sizes of a collection of 
several thousand clumps, in the fixed potential of the initial cloud and subject to 
the feedback processes described above. At the same time we model the coupled 
evolution of a spherically symmetric wind cavity and HII region. While we do not 
account for clump collisions, we note that their initial mass spectrum is probably 
created by a steady-state balance of collisional fragmentation and agglomeration. 
Our model does yet include a physical mechanism for star formation and so we 
investigate the response of our fiducial gas clouds to different imposed steady 
star formation rates, such that 50% of the cloud would be converted to stars after 
3 Myr. Stars are created from gas in the innermost clumps at these rates and 
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artificially added to the central cluster. We start our models by instantly forming 
~ 250 Mq of stars, equivalent to an ionizing luminosity of S49 = 1. We examine 
the time taken for the HII region to reach 90% of the initial cloud radius, at which 
point star formation is assumed to cease. Photoionization gradually destroys our 
model clouds, A, B and C, as shown in figure [f]. 

5 Discussion 

The turbulent, clumpy and self-gravitating nature of our clouds impedes their 
destruction and allows star formation to proceed to higher efficiencies. For ex- 
ample, an HII region in a uniform, quiescent cloud with the same mean density 
and initial star cluster as model A, would destroy the cloud in ~ 3 x 10 5 years, 
even with no additional star formation. However, with most of the cloud mass 
in dense clumps with velocities set by virial equilibrium, the process of "turbu- 
lent mass loading" confines the ionized gas for 1.5 - 2 Myr, even though by this 
time, five times the initial stellar mass has formed and feedback is correspond- 
ingly greater (model A). Reasonable models for clouds forming SSCs survive up 
to 3 Myr (models B and C). The fraction of the HII region occupied by the wind 
cavity increases with star cluster mass. 

After a few Myr our neglect of Wolf-Rayet winds and supernovae becomes 
important. Furthermore, our use of a fixed gravitational potential, our assump- 
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1 

time (Myr) 



Fig. 1. Star cluster formation and cloud destruction at constant star formation rates, 
<j>50- A: M c = 4000 M (Galactic case); B: M c = 4 x 10 4 M ; C: M c = 4 x 10 s M . 
Star formation efficiencies are quoted when Rst = 0.9i? c or after 3 Myr. 
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tion that all massive stars are at the cloud center and the absence of protostellar 
winds in our model cause us to overestimate the cloud destruction time. Never- 
theless, our results are qualitatively consistent with observed ages of the youngest 
optically visible SSCs, which are ~ 5Myr p2fl . Younger clusters are still embed- 
ded in dense gas 0,|l5)- R136 in 30 Doradus, which is considered a small SSC, 
has recently dispersed its gas and is ~ 1 — 2 Myr old |T^] . 

With steady star formation rates over the cloud lifetime, our results imply 
that star formation efficiency increases with initial cloud mass. The high efficien- 
cies apparent in our SSC models allow for the creation of bound clusters Jll| , 
even in the presence of vigorous feedback. Since the mass loss is gradual, efficien- 
cies as low as ~ 30% may result in loosely-bound clusters. Our Galactic model 
(A) is close to this limit, though the inclusion of additional feedback processes, 
such as protostellar winds, will reduce the efficiency. 

We plan to extend our models to include additional feedback processes, such 
as protostellar winds, Wolf-Rayet winds and supernovae, and a physical star for- 
mation mechanism (e.g., photoionization regulated Q) in the neutral clumps. 
We shall compare our models to observations of extra-galactic compact HII re- 
gions (e.g., 0). We hope to predict minimum cloud masses and densities re- 
quired to form bound stellar systems and probe in more detail the differences 
between Galactic and super star clusters. 
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